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By taking full advantage of the dynamic property imposed by the detailed balance condition, we 
derive a new refined unified fluctuation theorem (FT) for general stochastic thermodynamic systems. 

This FT involves the joint probability distribution functions of the final phase space point and a 
thermodynamic variable. Jarzynski equality, Crooks fluctuation theorem, and the FTs of heat as 
well as the trajectory entropy production can be regarded as special cases of this refined unified FT, 
and all of them are generalized to arbitrary initial distributions. We also find that the refined unified 
FT can easily reproduce the FTs for processes with the feedback control, due to its unconventional 
structure that separates the thermodynamic variable from the choices of initial distributions. Our 
result is heuristic for further understanding of the relations and distinctions between all kinds of 
FTs, and might be valuable for studying thermodynamic processes with information exchange. 
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I. INTRODUCTION 

There have been many great progresses in nonequi¬ 
librium statistical physics of small systems in the past 
two decades. Compared with classical statistical physics 
where relative thermal fluctuations are generally Gaus¬ 
sian and vanishingly small, fluctuations become much 
more prominent in small systems that undergo processes 
arbitrarily far from equilibrium due to some nonequilib¬ 
rium external drivings. Despite of the complexity that 
originates from the arbitrariness of the driving protocols, 
these fluctuations turn out to satisfy some strong, use¬ 
ful and elegant universal properties, exactly depicted by 
the fluctuation theorems (FTs) [1-6]. For example, the 
Jarzynski equality (JE) [1], which connects the work done 
in nonequilibrium processes and the free energy difference 
of the system at the initial and the final stages, is an inte¬ 
gral fluctuation theorem (IFT) of work, while its stronger 
version, the Crooks fluctuation theorem (CFT) [3], is a 
detailed fluctuation theorem (DFT) of work. It should 
be emphasized that both the JE and the CFT require 
the system initially prepared in a canonical distribution 
and the validity of the detailed balance (DB) condition. 
As a more universal IFT, the entropy production iden¬ 
tity (EPI) [5] holds for arbitrary initial distributions even 
without the DB condition. It is worth mentioning that 
a DFT of heat has been discovered quite recently, where 
the distribution of the initial state in the phase space is 
required to be uniform rather than canonical [6]. 

Here come several fundamental questions: Why are 
there strict requirements for the distributions of the ini¬ 
tial state in these FTs, and can they be released? Can 
these existing FTs be traced back to the same root? 
Besides the FTs of the work, the heat, and the trajec¬ 
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tory entropy production, are there any other FTs as¬ 
sociated with other thermodynamic variables? These 
questions were partially answered by the unified FTs ex¬ 
plored by Seifert [7, 8]. For a stochastic system under¬ 
going a nonequilibrium process within the time interval 
[0, r], driven by a temporally varying work parameter At 
(0 < t < r), the unified IFT reads [7, 8] 

(e-«) = = I , (1) 

To(ro) 


where the trajectory-dependent functional i?[rt] = 
ln[T’(rt)/:P(ft)] = As„, - lnpo(r]l) -b Inpo(ro)- Here 
Asm is the entropy production in the medium [5]; p“(r) is 
the auxiliary final distribution in the phase space; Po(F) 
(Po(r)) is the initial distribution of the system in the 
phase space for the forward (backward) process; ViTt) 
('P(ft)) is the probability density of a trajectory Ft (Ft) 
in the trajectory space due to the protocol At (time- 
reverse protocol At = Xr-t); ft = F[i)_t is the time- 
reverse trajectory of Ft with F^ to be the time-reversal 
of the phase space point F (e.g., for the underdamped 
Langevin dynamics. Ft = (rt,Pt)j then f| = (rt,—Pj) 
and ft = (rr-t, —Pr-t))- In fact, Eq. (1) is a unification 
of the IFT of heat, the JE and the EPI, since they can 
be respectively generated by setting the initial distribu¬ 
tion po(ro) and the auxiliary final distributions p“(r,-) 
to be both uniform, both canonical and the distribu¬ 
tions connected by the real dynamic evolution. For those 
processes the thermodynamic variable Sa with a definite 
time-reversal parity, say Sapt] = ect>S'Q[rt], Ca = ±1, we 
further have the following unified DFT [8] 


P{{Sa = egSg}) 

P{{Sa = Sa}) 


(e = So,}) , 


( 2 ) 


where P{{Sc,}){P{{Sa})) denotes the joint distribution 
of multiple variables Sa for the forward (backward) pro¬ 
cess [9, 10], and the right hand side (r. h. s.) means 
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the conditional expectation of e ^ when Sa = Sq- To 
generate the DFT of heat or the CFT, we simply choose 
the single odd parity quantity S to be the heat (Q) or 
the work {W), then set po(r) and po(r) to be both uni¬ 
form distributions or both canonical distributions. This 
choice leads to i? = jSQ or j3{W — AF) and thus makes 
the r. h. s. of Eq. (2) simply e~^‘^ or since 

R is uniquely determined by S and the conditional dis¬ 
tribution is a delta one. From this point of view, we 
can say that the FTs mentioned in the last paragraph 
do share the same root, which might be summarized as 
a combination of the microscopic reversibility (MR) and 
the dynamic property. The MR ensures the validity of 
Eq. (1) (and Eq. (2), though not obvious), since it is no 
more than the probability-normalization relation of all 
the time-reverse trajectories whose forward ones are of 
nonzero probability. In fact, the conventional ETs are 
inapplicable to absolutely irreversible processes because 
of the breaking down of the MR, as has been highlighted 
in recent investigations [11]. The dynamic property, such 
as the DB condition, is necessary to relate certain ther¬ 
modynamic variable to R by properly choosing the two 
initial distributions. This is important for endowing spe¬ 
cific physical meaning to the abstract identities (1) and 
(2) as merely the corollaries of the MR. 

On the other hand, despite of their universal valid¬ 
ity, the unified IFT (1) will be physically meaningless if 
R can’t be related to certain thermodynamic variables, 
while the r. h. s. of the unified DFT (2) is usually diffi¬ 
cult to either calculate or be given a transparent physical 
interpretation, unless R uniquely depends on S. Such en¬ 
tanglement between the thermodynamic variable and the 
initial distributions is the reason why the EPI holds for 
arbitrary distributions while the ETs of the heat and the 
work require specific initial distributions: the functional 
of entropy production naturally contains the distribution 
functions of the initial and final phase space points, but 
the heat and the work have nothing to do with them. 

If the initial distribution is not canonical (uniform) 
distribution, can we still construct a fluctuation theorem 
for the work (heat)? According to the above analysis, 
the answer seems to be no. However, in this article, we 
propose a new refined unified ET which achieves the sep¬ 
aration of the thermodynamic variable from the choices 
of the initial distributions. As a result of the DB con¬ 
dition, this refined unified FT is even more “detailed” 
than DFT (2), because it involves the joint distributions 
with the phase space point. In the studies of FTs, the 
first encounter of the joint distributions with the phase 
space point is in the Hummer-Szabo relation [2, 4, 12]. 
However, the Hummer-Szabo relation still requires the 
distribution of the initial state to be a canonical distri¬ 
bution, and is valid for the work only. In our current 
investigation, we extend the FTs of work and heats to 
an arbitrary initial distribution and even to other vari¬ 
ables, such as the entropy production. The cost is that 
we need to know the joint distribution function with the 
phase space point, such as Pt{W,T), which is more de¬ 


tailed than the usual distribution function, such as the 
work distribution Pt{W). We show that the new refined 
unified FT can generate many existing FTs as well as 
many new FTs that are not known previously by choos¬ 
ing proper initial states. It may also be of potential values 
in investigating information thermodynamics, where the 
initial distributions might be quite irregular due to the 
extra distribution rectification by the information. 

We notice that similar problems, i.e., the FTs for arbi¬ 
trary initial states, are investigated in a recent work [13]. 
However, the systems considered therein are discrete, au¬ 
tomatons (not externally driven) and may not satisfy the 
DB condition, and the FTs are associated with the local 
currents in a general dynamic network. Thus both as¬ 
pects of their study are complementary to ours. This 
issue is also discussed in Ref. [14], but is not its main 
focus. 

This paper is organized as follows. In Sec. H, we derive 
the refined unified FT by analyzing its close relation to 
the DB condition for general stochastic thermodynamic 
systems. In Sec. HI, we reproduce the existing FTs and 
generate some new ones of the work, the heat and the 
entropy production as some examples of the refined uni¬ 
fied FT. Some applications of the refined unified FT are 
explored in Sec. IV. In Sec. V, we summarize the paper. 
Six appendices are added at the end of the paper, which 
we think are helpful for understanding the details of the 
paper and the relevant issues. 


II. MAIN RESULT AND ITS DERIVATION 

A. Main result — the refined unified FT 

We focus on general stochastic thermodynamic sys¬ 
tems with the DB condition. The system is coupled to a 
heat bath with the inverse temperature /3, thus undergoes 
isothermal processes within the time interval [0, t] due to 
an external driving protocol denoted by At, t G [0, r]. We 
consider a trajectory-dependent thermodynamic quan¬ 
tity A [Ft] in the following form 

A[rt] = /3g[rt] + a,(ro-ao(ro). (3) 

Here the heat functional g[rt] = — dtTtdrUtiTt) 
{Q > 0 corresponds to the release of the heat from the 
system to the heat bath), with t/f(r) to be the energy 
of phase space point F at time t; at(r) can be an arbi¬ 
trary time-dependent function of the phase space point 
F, for which we can generally define its time-reversal 
at(r) = aT-t(rl) [15] and correspondingly 

Apt] =/3Q[rt]+a,(r,)-ao(ro) . (4) 

All over the paper, if not particularly indicated, we al¬ 
ways stipulate that the energy possesses the property of 
the time-reversal invariance, i.e., t/t(rl) = C/t(r). Our 
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main result reads 

/g^dm(-^,r)po(rt)e-(n 

/g^drp,(A,r)po(rt)e“^(r) ® 

where the accessible phase space at time t {t — t) for the 
forward (backward) process is denoted by &t, thus the 
integral of F in the denominator (numerator) goes over 
the whole accessible phase space of Ft- (Fq or Ft- = fJ) 
[16]; Po(r) (po(r)) denotes the distribution of the initial 
state in the phase space for the forward (backward) pro¬ 
cess; Pt{A, F) F)) is the joint distribution function 

of the thermodynamic variable A (A) accumulated until 
time T and the final (at time r) phase space point, start¬ 
ing from the initial distribution po(r) (po(r)) and driven 
under the protocol A* (At = At— t)- Eq. (5) is valid for 
arbitrary initial distributions Po(r) and po(r), and arbi¬ 
trary state functions at(F). Also, similar to the JE and 
the GET, it is valid for an arbitrary protocol At and an 
arbitrary driving time r. 

B. Dynamic property equivalent to the DB 
condition 

In order to demonstrate the close relation between our 
main result (5) and the DB condition transparently, we 
use a somehow complicated method to carry out the 
derivation, though a relatively simple but much more 
mathematical proof is available by using the path integral 
approach (see Appendix A). 

We first write down the Fokker-Planck equation (FPE) 
of the stochastic system 

diPtiX) = Apt(F) , (6) 

with the generator Ct to be a time-dependent linear 
operator corresponding to the protocol At and only 
acting on F. For later use, we define the trans¬ 
pose operator of Ct (denoted by Cf), which satisfies 
/e* ^r 5 (F)£t/(F) = J^^drf{T)Cfg{r) for arbitrary 
normalizable (i.e., | f dF/(F)| < -|-oo) functions /(F) and 
( 7 (F) defined in the phase space. One can see that the 
normalization of Pt (F) will impose the property C^ 1 = 0 
to Ct if we integrate the FPE ( 6 ) over F, though this 
is not rigorous for the systems with infinite phase space 
where 1 is not normalizable. We can also define the time- 
reversed operator of Ct (denoted by C\). It is obtained 
by adding minus signs to all the components in Ct with 
the odd time-reversal parity (e.g., p ——p, dp —dp, 
p is momentum). 

In terms of the generator, the DB condition manifests 
itself in the following algebraic symmetry 

^ jrT ( 7 ) 

Such dynamic property comes directly from the com¬ 
bination of the original definition of the DB condition 
g-/3£/t(ri)^^^P^ _>• F 2 ) = e~^^*‘'^^'>Wt{T\ —>■ f|) (for any 


Fi,F 2 € &t) and the transition rate formula Wt(Fi —>• 
r 2 ) = /g^ dr^(F — F 2 )£t( 5 (F — Fi) [17]. Conversely, the 
original DB condition follows if Eq. (7) is previously as¬ 
sumed, thus the dynamic property is actually equivalent 
to the DB condition. We emphasize that Eq. (7) should 
be valid at any time t € [0 ,t] despite the work parame¬ 
ter is temporally varied. This is the straightforward dy¬ 
namic generalization of the common static version (where 
the work parameter is fixed), which can be found in the 
standard literature [17]. We also want to mention that 
Eq. (7) ensures when we suddenly stop the external driv¬ 
ing at time t, the system will always relax to the canon¬ 
ical distribution determined by At owing to the property 
£^1 = 0. In addition, Eq. (7) is much more stronger than 
merely imposing = 0 , that’s why the balance 

is called detailed. 

As an example, it can be checked that Eq. (7) is true 
for the Langevin dynamics (in both the underdamped 
and the overdamped regimes, see Appendix B). It is also 
notable that a discrete but more general (may be without 
the DB condition) version of Eq. (7) has appeared in 
Ref. [18] (see Eq. (3.24) therein). 

C. Derivation of Eq. (5) based on the characteristic 
function 

According to the FPE ( 6 ) as well as the definition of 
the functional A (3), it can be proved (see Appendix C 
for a more general formula) that the joint distribution 
function Pt{A,r) satisfies the following equation of mo¬ 
tion 

dtPt{A,T) = je-[“‘(r)-/3(7t(r)]a^£^g[adr)-/3(7dr)]a^ 

- dtat{r)dA}PtiA,r) . 

Notice that the operator acting on Pt{A, F) on the r. h. s. 
only contains dA but is independent of A. It is natural to 
(partly) perform an integral transformation, as is always 
done when we construct the Feynman-Kac formula (e.g., 
at(F) = dUt{T) in Eq. ( 8 )). We take the inverse Fourier 
transformation to change Pt{A,r) into its characteristic 
function Gt( 7 i,F) = dAe*'^^Pt(A, F). For further 
simplicity, we define the modified characteristic function 
G^{p,T) = . Based on Eq. ( 8 ), it is 

found that G]"( 7 x,F) satisfies the following equation of 
motion 

dtGTip,r) = Ctip)GTip,r), (9) 

where Ct{p) = For the backward 

process, we define Ct{p) = £r-t(/i), so the equation 
of motion of the modified characteristic function of the 
backward process, G["( 7 i,F) = F), can be 

expressed as 

dtGT{p,T)=Ct{p)GT{p,T). (10) 

In terms of Ct{p), the relation (7) can be rewritten as 

^Kd + i) = ^T-ti-p) ( 11 ) 
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by using the identity £T^-it^pu,{r)_ 

order to make full use of Eq. (11), we write in the fol¬ 
lowing special forms the equation of motions for both the 
forward (9) and the backward (10) processes 

+ z,rt) = 4+ ,)G--(^ + z, rt), (12) 


- r) = Cr-t{-^i)G^.t{-^^, r). (13) 

Combining these two equations of motion with Eq. (11), 
after a straight forward calculation, we deduce that 


dt [ drGr(M + brt)G™ t(-A*,r) = o. (i4) 
Jet 


This means the integral dTGY^{^ + i,V^)G'^_^{^,T) 
is a conserved quantity during the dynamic evolution. 
Particularly, at the terminal time points [t = 0, r), we 
have 


[ drG™(Az + *,r)G^(-M,rt) 

Je^ 

= [ drG™(/z + z,rt)G™(-Az,r). 
Jen 


(15) 


After substituting the expressions of the modified char¬ 
acteristic functions into Eq. (15), we obtain 


[ drG,(/r + z,r)po(r^)e“-(i’) 
Je-r 

= [ drG,(-/z,r)po(r^)e“-(r) ^ 

Jeo 


(16) 


By taking the Eourier transformation of Eq. (16), we fi¬ 
nally obtain our main result (5). 

It is worth mentioning that a similar method based on 
the symmetry of the generator has appeared in a recent 
work [19] (see Eq. (11) here and Eq. (10) therein). But 
their focus was on the Bochkov-Kuzovlev equality [20] for 
the open quantum systems described by Lindblad mas¬ 
ter equations. Therefore, we believe there is a quantum 
generalization of the main result (5) (and its extensive 
corollaries) for at least Lindblad-type open quantum sys¬ 
tems, which we leave for our future work. 


III. REFINEMENT OF THE EXISTING FTS 
AND GENERATING OF NEW FTS 

A. refined FT of work 


The simplest specialization of the refined unified FTs 
(5) is when we choose at(r) = j3Ut{T). In this manner, 
both the quantity A (3) and A (4) turn out to be the 
dimensionless work fJW, due to the first law of thermo¬ 
dynamics at the level of individual trajectories [25] . The 
refined FT of work reads 


j^jrPri-w,r)po{r^)eP^o(^'^ 

J^JTPriW,r)Poir^)edUo(ri) 


Similar to the GET, this relation is valid for an arbi¬ 
trary driving protocol At and an arbitrary driving time 
T. What is more, this relation is more general than the 
GET because it is valid for arbitrary initial distributions 
Po(r) and po(r). Obviously, if we want to construct the 
existing work FTs from Eq. (17), the choices of the dis¬ 
tributions of the initial state in the phase space should be 
the canonical ones for both the forward and the backward 
processes. That is, _po(r) = = e“^^“(^V-^o(/3) 

and po(r) = Po'^(r) = where the par¬ 
tition functions ^o(/3) = and Zo(/3) = 

/g For such choices, Eq. (17) reduces to 

the well-known GET [3] 


P^(-W) ^ 

Pr{W) 


(18) 


where Pr{W) = J^^dTPriW^r) and Pr{W) = 
/g^ dTPr{W,T), AF = -/3-1 ln[Z^(/3)/Zo(^)] is the free 
energy difference [29]. The integral version of the GET 
(18) is the celebrated JE [1] 

(e-'^(^-A^)) = 1 . (19) 


We can also easily reproduce the Hummer-Szabo relation 

[ 12 ] 

„-/3c/^(r') 

(.(r.-rOe-^-) = ^^, (20) 

as long as we set the initial distribution for the forward 
process to be a canonical distribution pQ®(r) and the ini¬ 
tial distribution for the backward process to be a (5 dis¬ 
tribution po(r) = <^(r — E^^) in Eq. (17). 

The discrete version of Eq. (17) is 

Sm Pr{-W, m)po(z7it)e^^E ^ 

EnPr{W,n)poinr)e^^- e , ^ j 

where E" {E'P) is the energy of the state n {fh) at time t 
for the forward (backward) process; denotes the time- 
reversed state of n, e.g., the spin-down state if n denotes 
spin-up. Gertainly, this relation holds for discrete-level 
stochastic systems. We would like to emphasize that 
Eq. (21) is also valid for an isolated quantum system (not 
necessarily with time-reversal symmetry), where the ini¬ 
tial density matrix is generally po = Sm 7 'o(nz)|'nz)(m| 
(and ffo = X]nPo(zz)|h)(n| for the time-reversed process), 
as long as we admit the two-point projection measure¬ 
ment definition of quantum work (see Appendix D). Fur¬ 
thermore, (3 in Eq. (21) can be arbitrarily chosen and may 
even be a complex number, since the quantum system is 
isolated from any heat bath. 


B. generating new FTs of work 

In addition, from the refined FT of work (17) one 
can derive several FTs of work that were previously not 
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known to researchers in this field. For that purpose, let us 
first choose both po(r) = i5(r—Fg) andpo(r) = i5(r—F'l) 
to be 5 distributions. So that Eq. (17) becomes 




pT(r't) 

Po^(ro) 


Pr{-w,Tl\r'^) , 
( 22 ) 


where PT(kF,F'IF q) is the conditional joint probability 
distribution. It means the sum of the probabilities of 
all the forward trajectories that end at F', and the work 
accumulated along each of these trajectories are equal to 
W conditioned on the given initial phase space point Fg. 
P.,.(—IE, rJlF'l) can be understood in a similar way but 
for the backward process. This relation can be regarded 
as a generalization of CFT to initial 6 distributions, and 
has previously been obtained in Refs. [21, 24]. We do 
integral with respect to W on both sides of Eq. (22) and 
obtain 


(<5(r. - F') 


g-/3(lV-AF) 



pT(r't) 

Po^(ro) 


Pr(rJir't) 


(23) 


where (E[Ft]|Fo) indicates the ensemble average of the 
functional E[rt] over all trajectories that start from the 
phase space point Fg, and Pr(Fj|Fh) is the conditional 
probability distribution in the phase space for the back¬ 
ward process. It describes the final probability distribu¬ 
tion of the backward process at fJ given the initial state 
at F'l. This relation (23) can be regarded as a general¬ 
ization of the Hummer-Szabo relation (20) to the S initial 
distribution. If we further do integral with respect to F' 
on both sides of Eq. (23), we obtain a JE-like FT 


^g-/3(W-AF) 



PT(rg) 
Po’(ro) ’ 


(24) 


where Pr(r) is the final probability distribution in the 
phase space evolved from the initial canonical distribu¬ 
tion pg'^(r) in the time-reversed process. This relation 
can be regarded as the generalization of JE to the <5 
initial distribution. Similar mathematical relations have 
been obtained in Ref. [22] for isolated classical systems 
and in Ref. [23] for open quantum systems. But in both 
these two references, their focus are on different prob¬ 
lems, and the mathematical relations are not interpreted 
in this way. 

Besides the initial 5 distribution, the JE (19) and the 
Hummer-Szabo relation (20) can actually be extended to 
an arbitrary initial distribution in the forward process. 
We multiply an arbitrary initial distribution pg(rg) to 
both sides of Eq. (22) before we do the integral with re¬ 
spect to Fg, and we obtain the extended Hummer-Szabo 
relation for an arbitrary initial distribution pg(F) 


(^(5(F^ 

= [ dFp,(rt|F't)p[;«(r't) 

Jen 


poT) 

po(n 

pIHt) 


where (•■•)pp(r) represents the average over all trajecto¬ 
ries when the initial distribution is Po(r) for the forward 
process. If we further do integral with respect to F' on 
both sides of Eq. (25) we obtain the generalized JE for 
an arbitrary initial distribution 


-I 3 {W-AF) 




(26) 


In this equality, the initial distribution of the forward 
process can be an arbitrary distribution pg(r), while the 
initial distribution of the backward process must be a 
canonical distribution. It is worth mentioning that in 
the FT of the total entropy production along individual 
trajectories [5], the initial distribution can be an arbi¬ 
trary distribution (38). It was assumed previously that 
only when the initial distribution is a canonical distribu¬ 
tion (globally thermal equilibrium) or a partially thermal 
equilibrium distribution [26, 27] can one construct the FT 
of work. Here we generalize the FT of work (mainly JE 
(19), CFT (18) and Hummer-Szabo relation (20)) to an 
arbitrary initial state. The generalized JE for an arbi¬ 
trary initial state (26) may have potential applications 
in free energy recovering via nonequilibrium work mea¬ 
surement and numerical free energy computation. 


C. refined FT of heat 


Another specialization of the refined unified FT (5) is 
the FT of heat when at(F) equals to a constant (e.g., 
zero), thus A = A = (3Q. The refined FT of heat reads 

Je„JrF(-o,r)po(rt) ^ 
S^dTP,{Q,r)p„{Tt) ■ ' ' 

Similar to the DFT of heat [6] , this relation is valid for an 
arbitrary driving protocol At and driving time t. More¬ 
over, this relation is more general than the DFT of heat, 
because it is valid for arbitrary initial distributions Pg(r) 
and pg(r). If we want to reduce the joint distribution 
function to the heat distribution function, the only choice 
is to make both pg(F) and pg(r) independent of the phase 
space point, i.e., the uniform distribution, as was found 
in Ref. [6]. However, a uniform distribution can never be 
truly realized in continuous systems, where the entropy 
has no upper bound due to the infinite volume of the 
phase space. So we have to consider finite-level systems, 
such as spin systems. Similar to Eq. (21) , the discrete 
version of Eq. (27) is 


Srn-Pr(-Q,m)pg(TO^) ^ pQ 

Y.uPr{Q,n)po{n^) 


Suppose that the system has totally N states. By setting 
pg(n'l') = po{m') = 1/N ( maximum entropy state) in 
Eq. (28), we obtain 


Q) ^ - 0 Q 

PriQ) 


( 25 ) 


( 29 ) 
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where Pr{Q) =Y.nPr{Q:n) and Pr{Q) =Y.rnPr{Q,rn) 
are respectively the heat distribution functions for the 
forward and the backward processes. Eq. (29) is nothing 
but the DFT of heat, whose integral version reads [6] 

(e-PQ) = 1 . (30) 

It is worth mentioning that Eqs. (29) and (30) are likely 
to be experimentally tested in analogy to the verification 
of the EPI for a finite-level system [28]. More generally, 
we can test Eq. (27) in whatever systems with finite phase 
space. 

One may imagine that the uniform initial distribution 
of the DFT or the IFT of heat might be approached via 
the limit /3' —>• 0 (T' oo), with fi' to be the inverse 
temperature of the initial canonical distribution. Hence, 
one may ask whether 


to > ,„4M_ 

/3'—>0 f5Q Pt{—Q) 


= 1 


(31) 


holds true even for a continuons system with its initial 
distribution to be po(r) = e~^ /Zo{P'). In fact, it 

has been demonstrated [6] that for a driven Brownian 
harmonic oscillator, if we rewrite Q in wrp , with w to be 
a positive quantity that characterizes the driving speed, 
and take the limit r —>■ oo (so that Q ^ oo for finite p) be¬ 
fore /3' —0, Eq. (31) will be indeed true. It might be an 
intrigning but involved snbject to investigate the valid¬ 
ity range of Eq. (31) for more general kinds of stochastic 
systems with infinite phase space. 


D. generating new FTs of heat 

Similar to JE for an arbitrary initial distribution, we 
can extend the integral FT of heat (30) to an arbitrary 
initial state. For the sake of a well-defined uniform distri¬ 
bution, we write down the discrete version for an iV-state 
system 


. (32) 

where Po(m) = 1/JV is the initial uniform distribution 
of the forward process, and the final distribution of the 
microscopic state for the backward process Pr{'m)) = 
f dQPr{Q,m'^) in Eq. (32) must correspond to the uni¬ 
form initial distribution. One can see that the r. h. s. of 
Eq. (32) usnally deviates from 1, if the initial distribu¬ 
tion of the forward process is not prepared in the uniform 
distribution. Eq. (32) is valid for an arbitrary protocol 
and an arbitrary driving time. What is more it does not 
require the initial distribution of the forward process to 
be the uniform distribution. Hence, if Eq. (27) can be re¬ 
garded as the generalizations of the detailed version of the 
heat FT [6] to an arbitrary initial distribution, Eq. (32) 
can be regarded as the generalization of the integral ver¬ 
sion of heat FT [6] to an arbitrary initial distribution. 


For the heat FT, one can also derive a Hummer-Szabo 
like relation. The initial distribution for the forward pro¬ 
cess and the reverse process are chosen to be the uniform 
distribution and the <5 distribution respectively. Let rrir 
be the final (at time r) state of the forward process, then 
this relation reads 

^ . (33) 

Furthermore, we can also extend this Hummer-Szabo like 
relation for heat (33) to an arbitrary initial distribution 

Po{m) 


{Sm^ne = X! ^^0 ■ (34) 

m 

Up to now, we have successfully generalized the FTs 
of work and heat to arbitrary initial distributions. For 
a comparison between previous results and our current 
results, please see Table I. 


E. Unified IFT and EPI 


Besides reproducing the JE, the GET and the heat 
FTs, the refined unified FT (5) can reproduce the EPI. 
Let’s again have a look at the main resnlt (5). We find 
that if aT(r) = — Inpo(rl) and aT-(r) = — Inpo(r^), 
which are sufficient to determine ao(r) = ar(r^) = 
— Inpo(r) and ao(r) = aT(r^) = — Inpo(r) (thus H.[rt] 
and ^[ft] are completely determined) dne to the defini¬ 
tion of the time-reversal at(r) = a,—t(r’^), Eq. (5) will 
simply reduce to 


PA-A) 

Pr{A) 


(35) 


Here Pr{A) = /g^drP^(H,r) and Pr{A) = 
J^^drPr{A,T). The quantity H[rt] = /3Q[Tt] - 
lnpo(r]l) + Inpo(ro) depends on the initial distributions 
of both the forward and the backward processes, which 
is qnite different from the cases of the heat and the work. 
The integral version of Eq. (35) reads 


(e 




PA^r) 

Po (bo) 


e-W) = 1 , 


(36) 


where the auxiliary final distribntion in the phase space 
p“(r) = po(r^) can be an arbitrary one, since the choice 
of Po(r) has no restriction. The generalized version of 
Eq. (36) to the cases without the DB condition, which is 
merely to replace f3Q with As^, is the unified IFT (1) 
we mentioned in the beginning. 

Despite of the fact that the quantity A defined in such 
way always satisfies the DFT (35), usually it cannot be 
related to any thermodynamic variable we are familiar 
with. To see this, we consider the average of A over all 
the trajectories 

(A) = (3{Q) + D[pr{r)\\po{TA] + {s)r - {s)o , (37) 
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where i:)[pi(r)||p 2 (r)] = / drpi(r) ln[pi(r)/p 2 (r)] is the 
Kullback-Leibler divergence between two probability dis¬ 
tributions _Pi(r) and _P 2 (r); {s)t is the ensemble-average 
entropy of the system at time t; PT(r) is the final dis¬ 
tribution in the phase space determined by the real dy¬ 
namic evolution. A special case in which A has a trans¬ 
parent physical meaning is when po(r^) = pT-(r). In 
this case the second term in Eq. (37) vanishes. For 
such a choice of the initial distribution for the back¬ 
ward process, A is the trajectory-dependent total en¬ 
tropy production Astot for the forward process, and 
(A) is in consistency with the ensemble-average value 
{Astot)- However, since dr(r) has been determined by 
po(r) as — Inpo(rf), po(rf) is usually different from the 
final phase space point distribution Pr(r) starting from 
the initial distribution po(r) = Pr(rf) and driven by the 
time-reversed protocol [30]. As a result, the functional 
A[ft] = /3Q[rt] —lnpo(ff )-|-lnpo(fo) cannot be regarded 
as the total entropy production of the time-reversed tra¬ 
jectory Ft. Hence, even when we choose po(rf) = Pri^), 
Eq. (35) cannot be regarded as a detailed EPI. On the 
other hand, Pt{A) is still a normalized distribution func¬ 
tion. So we can obtain the EPI [5] by setting A to be 
Astot without any problem 

(e-^"“>‘) = I . (38) 

Nevertheless, what if Po(r^) happens to be the final 
distribution of the time-reversed process? In fact, such 
kind of specific distribution, named by echo state, has 
been recently studied by Van den Broeck and collabora¬ 
tors [31], which clarifies the initial conditions of both the 
detailed EPI and the other DFTs in Ref. [32]. Similar to 
the discrete case that Van den Broeck et al discussed, for 
any given FPE (6), the corresponding echo state 
for the detailed EPI can be uniquely determined by solv¬ 
ing 

p“^°(r) = T{eio }r{e^o ‘^‘^*}p“^°(r), (39) 

where T{...} denotes the time-ordered product. For such 
choice of the initial distribution, the detailed EPI 


Pt{ Astot) ^ g-Astot 

Pr{Astot) 


(40) 


holds unambiguously. This can be regarded as a gener¬ 
alization of the detailed EPI for steady state [5], where 
Ct should be time-independent, i.e., the work parameter 
should be fixed. 


IV. APPLICATIONS OF THE REFINED 
UNIFIED FT 

A. Calculating distribution functions of work and 
heat for arbitrary initial states 


Po(r) = e “oHl/Jg dTe obtain 

P^{A) = e^ [ / driV(-A,r)po(r^')e“"(^^ 

J ©T J Go 

(41) 

This equation implies that Pr (A) can be easily computed 
for arbitrary initial distribution po(r) as long as we know 
the joint distribution function Pr{A^T) corresponding 
to the particular initial distribution Po(r) for the time- 
reversed process. Specially, if Po(r) is a delta function, 
Eq. (41) will become 

P^(Alro) = e^+“'>(^°)p^(-A,rJ) / (42) 

Je-r 

where Pr(Ajro) denotes the distribution of A in condition 
that the initial state is known to be Pq. In fact, Eq. (41) 
and Eq. (42) are equivalent to each other, similar to the 
equivalence between the dynamic property (7) and the 
DB condition, due to the additivity of the probability 
for exclusive events, Pt{A) = c?rPr(^|r)po(r). As 
two examples, we can calculate the work and the heat 
distributions for any initial distributions Po(r) {po(rn)) 
via the following two formulas 

P,{W) = [ dTPr{-W,r^)^^ , (43) 

Jeo PoA) 

f..((3) = e«y:P.(-<3,mt)^. (44) 

m -^0 V ) 

These relations provide an alternative approach to obtain 
the work (heat) statistics besides, e.g., directly solving 
the Feynman-Kac formula for the forward process with 
an arbitrary initial condition. Such approach may be 
advantageous in certain cases due to the specific initial 
condition of the backward process. 

As an example, let’s consider an overdamped breathing 
Brownian harmonic oscillator in one dimension, of which 
the generator in the FPE reads 

Pt = -d^{ktx + j3~^d^) , (45) 

7 

with the protocol chosen to be kt = + nt. According 

to Ref. [33], using the Gaussian ansatz, the Feynman- 
Kac formula (partial differential equation) of this model 
can be reduced to a set of first-order ordinary differen¬ 
tial equations, which can be easily solved numerically. 
Furthermore, if the initial state is the equilibrium state 
and the driving speed is slow, we can even perturbatively 
work out some analytical results, such as the mean work 
correction in the linear response regime. However, the di¬ 
rect perturbative analysis breaks down if the initial state 
deviates significantly from the equilibrium one. To by¬ 
pass the difficulty, we make use the following equivalent 
form of Eq. (43) 


Let’s come back to the main result (5). If we only 
set the initial distribution of the backward process to be 


G,{p) = 




Po(n 

Po"(r) 


(46) 






where G'T-(/r,r) has an exact perturbative solution ow¬ 
ing to its equilibrium initial state. In particular, if we 
choose po{x) = (/3'fco/27r)2e“^ the canonical dis¬ 

tribution of another temperature using the conclu¬ 
sions in Ref. [33], we will obtain a Gaussian distribution 
in the slow driving limit, of which the mean reads (see 
Appendix E for details) 

while the variance cr^ is the same as that for the equi¬ 
librium initial state. Here the free energy difference 
AF = In y^kr/ko and kr = kg + kt . A notable inference 
follows that for an expanding process {kr > kg), {W) will 
be less than AF when -|- l)/3“^/2, which 

implies the necessity of an equilibrium initial state (with 
the same temperature of the heat bath) for the validity of 
the maximum work principle. Actually, the above anal¬ 
ysis can be generalized to arbitrary slowly driven over¬ 
damped Langevin systems, whose Feynman-Kac formu¬ 
las have a unified analytical treatment [34]. 

B. Deriving FTs for feedback control processes 

Although the canonical distribution is much more com¬ 
mon than other ones in real cases, the initial distribution 
corresponding to a specific protocol can be rather irreg¬ 
ular in feedback control processes. In the extreme cases, 
we can exactly measure the initial phase space point and 
then choose its unique protocol. More generally, we as¬ 
sign a protocol to a region in the phase space, while the 
region can be arbitrary in principle. According to these 
observations, we infer that our main result may have po¬ 
tential applications in thermodynamics with the feedback 
control. 

Actually, we can easily derive one of the Sagawa-Ueda 
equalities [35] based on our main result or its side prod¬ 
ucts. We apply Eq. (43) to a protocol \\ which corre¬ 
sponds to the measurement result y 

Pr{W\y) = [ dTPri-W,T\y)P^ff^ . 

Jeo PoP^) 

(48) 

Here the conditional initial distribution in the phase 
space po(r| 2 /) satisfies 

po{^\y)p{y) = p( 2 /|r)Po'^(r), (49) 

as long as the system is initially prepared in a thermal 
equilibrium state, with p{y) = dTp(y\T)pQ‘^(T) to be 
the probability density that the measurement result turns 
out to be y. p(y|r) is the conditional probability density 
that the measurement result is y conditioned on the phase 
space point E. For simplicity, we assume that the mea¬ 
surement has the property ^(j/'flrl) = p( 2 /lr) [35]. This 
relation is obviously satisfied for the error-free measure¬ 
ment p{y\T) = ^(j/ —r), and the Gaussian-error measure¬ 
ment p(ylr) = (27r(T^)“5e“l^“'"l , with E = (r,p). 


Combining Eq. (48) and Eq. (49), we obtain 

= f dTPri-W,T\y)piy\T^) . 

JGo 

(50) 

After integrating Eq. (50) over both y and W, we get its 
integral version 

{e-P(W-AF)^ = y , (51) 

where 

77 =/ dydrp{y'‘\T)pr(T\y) = f dyp{y^\y) . (52) 

J Go J Go 

Here Pr{T\y) = f dWPr{W,T\y) is the final distribu¬ 
tion of the phase space point for the backward process, 
and p{y'\y) is the probability density that the final phase 
space point is measured to be y', after being driven by 
the protocol X\ from the canonical time-reversed initial 
distribution. 

The other main result in Ref. [35] reads 

{e-P(w-AF)-i^ = 1 , ( 53 ) 

where I is the initial mutual information between the 
system and the measurement device. However, in both 
Eq. (51) and Eq. (53) it is required that the initial distri¬ 
bution must be a canonical distribution. To construct a 
more general information-involved FT, Sagawa and Ueda 
proposed 

^g-<T+A7) ^ ^ ( 54 ) 

in Ref. [36] as a generalization of the EPI. Here AI is the 
difference of the initial and the final mutual information; 
a = PQ + As and As = — InpT-(ri-) + lnpo(ro)i with 
Pt-(E) = J dypr{T\y)p{y). The above two Sagawa-Ueda 
equalities (53) and (54) are actually contained in the uni¬ 
fied IFT (1), if we use the point of view in Ref. [36] to 
regard the combination of the original system and the de¬ 
vice as a composite system (see Appendix F). However, 
the counterpart of Eq. (51) as another generalization of 
the EPI seems to be unexplored so far. Now we can de¬ 
rive it quite straightforwardly in analogy to the deriva¬ 
tion of Eq. (51). Following the same procedure as that 
from Eq. (48) to Eq. (50), we obtain 

p{y)Pricr\y)e~'" = [ drPri-(7,T\y)p{y\T'^) , (55) 

Jeo 

which indicates 

(e--^) = 77 . (56) 

The expression of 77 is exactly the same as Eq. (52), but 
p{y'\y) must be associated with the initial distribution of 
the time-reversed process, which is the real time-reversal 
of the forward final distribution, i.e., _Po(r) = Pr(r^)! 
instead of a canonical one. 
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TABLE I: Comparison between previous results and our results on the requirements on the initial state for different FTs. Please 
note that in previous results only distribution functions of certain thermodynamic variables, such as Pt{W) and Pr{Q), are 
required. But in our current results, more detailed joint distribution functions, such as Pt{W,V) and PT(Q,r), are required. 


Fluctuation Theorems 

Previous requirements on the 
distributions of initial state 

Our requirements on the 
distributions of initial state 

Work FTs 

CFT [3] 

JE [1] 

Hummer-Szabo relation [12] 

canonical distribution (18) 
canonical distribution (19) 
canonical distribution (20) 

arbitrary distribution (17) 
arbitrary distribution (26) 
arbitrary distribution (25) 

Heat FTs 

Detailed heat FT [6] 
Integral heat FT [6] 
Hummer-Szabo like relation 

uniform distribution (29) 
uniform distribution (30) 

N/A 

arbitrary distribution (27) 
arbitrary distribution (32) 
arbitrary distribution (34) 

FTs with feedback control 

Sagawa-Ueda equalities [35, 36] 

canonical distribution (51),(53) 
arbitrary distribution (54) 

arbitrary distribution (56) 


As an illustrative example of Eq. (56), we consider 
the Szilard engine [37], as was used in Ref. [35] to il¬ 
lustrate Eqs. (51) and (53). If the process is reversible 
and the measurement is error-free, we will always have 
As = 0 and cr = /3Q = —In 2 (delta distribution), so 
that {e~'^) = 2. On the other hand, if y = ^ (r), i.e., the 
particle is found in the left (right) half chamber, \\ will 
be the rightward (leftward) expansion. Thus will be 
the leftward (rightward) compression. It is easy to see 
that p(^l^) or p(rjr) (^f = ;, rf = r, since ^ or r is related 
to the position that is invariant under time-reversal) is 
simply 1 , owing to the fact that the particle must be al¬ 
ways in the left (right) half chamber after the leftward 
(rightward) compression. Hence, 77 = p{l\l) + p{r\r) = 2, 
and Eq. (56) is indeed valid for such feedback control 
processes. What’s more, even if the process is far from 
the quasistatic one, we will still have rj = 2 due to the 
same analysis. This result generally implies — (cr) < In 2, 
which is one aspect of the Landauer’s principle [38, 39]. 


V. CONCLUSION 

In summary, we propose a refined unified FT (5) under 
the DB condition. The refined unified FT is applicable 
to several thermodynamic variables, such as the heat, 
the work, the trajectory entropy production, and is even 
more refined than the DFT. Compared with the previ¬ 
ous unified IFT and DFT [7, 8 ], our refined unified FT 
(5) eliminates the entanglement between the thermody¬ 
namic variable and the choice of the initial distributions 
in the phase space, thus is physically more natural and 
comprehensible. In particular, our refined unified FT 
generalizes the FTs of the work and the heat, such as the 
JE and the CFT, to arbitrary initial distributions (see 
Table I), as well as generates several new FTs that were 
previously not known to researchers in this field, for ex¬ 
ample, the Jarzynski equality to an arbitrary initial state 
(26). We also revisits the validity of the DFT of entropy 
production [30-32] . The price that we pay for the gener¬ 
alization to arbitrary initial distributions is that we need 


to know the joint distribution functions with the phase 
space point. Due to such kind of generalizations, our re¬ 
sults might be valuable for the free energy recovering ex¬ 
periment and free energy computing as well as studying 
thermodynamics with information exchange, where the 
initial distributions are, to some extent, irregular. Based 
on the refined unified FT, we reproduce one (51) of the 
Sagawa-Ueda equalities and derive a new generalized EPI 
(56) for the feedback control processes. 
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Appendix A: Path Integral Derivation of the Main 
Result (5) 

In the path integral representation, the DB condition 
is reflected by the following relation 

1 1 ^(rtlTo) 

which is a crucial building block in the later derivation of 
the main result. Here P(rtjro) (■p(ftjfo)) is the condi¬ 
tional probability density of the trajectory Fj (F^) when 
the initial phase space point is known to be Fq (Fq) for 
the forward (backward) process. While Eq. (Al) has ap¬ 
peared in various references, such as its original discrete 
version in Ref. [40] , to make the paper self-contained, we 
briefly derive it from the basic trajectory definition of 
heat and the dynamic property (7). To do this, we first 
consider a sufficiently short time interval [t, t + dt] in the 
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forward process, during which 

p(rt+dt|rt) = [ drsir-rt+dt)ii+Ct+iidt)Sir-rt). 

(A2) 

Herep(rt,|rtj = dr5(r-rtjr{e^n 

is the conditional probability density that the phase space 
point at t 2 is starting from Fj^ (a delta distribution) 
at ti(< ^ 2 ), and the terms of the order of magnitude of 
(dt)^ are neglected. Similarly, for a short time interval 
[t — t — dt,T — t] in the backward process, we have 


Then we use the definition of A (3) to replace /3Q[Ft] + 
ariXr) with A[Ft] + ao(Fo), which leads to 


(A9) 


xpo(ro)d(A[Ft] - A) . 

Again we use the property of the delta function, getting 

^P[Fi]e“'>(i’«)iP(ft|fo):Po(rt) 
xpo(ro)<5(A[Ft] - A) . 


L = e^ 


(AlO) 


p{Tr-t\Tr-t-dt) 




dTd{T - f,_t)(l + C,_t_fdt)S{T - Tr-t-dt) 
dT6{T - Ft)(l + Cl^^dt)6{T - Ft+dt) . 

(A3) 


According to the definition of the time-reversal, it can 
be checked that A[ft] = —A[Ft], e.g., Q[ft] = —(5[Ft], so 
that 


L = e^ [ V[rt] [ dFe“-(i’-)iP(ft|fo):Po(fo) 

J dso (All) 

xpo{rl)S{A[rt] + A)6(Tr -T) . 


Combining Eqs. (A2) and (A3) with the dynamic prop¬ 
erty (7), we deduce that 

1^ MTt+dtlr^ ^ _pg jj (F,, .t)f,, dtdf . 
PiTr-tlTr-t-dt) 

(A4) 

Due to the Markovianity of the dynamics, the probability 
density of a trajectory can be expressed as P[Ft|Fo] = 
lim/f^oo rifci Therefore, we obtain 

T’[Ft|Fo] p(F K-k+i^\T 

= -p [ dtdrU{rt)tt = PQiVt] 

Jo 

Now we start the actual derivation of Eq. (5). Re¬ 
calling the definition of Pt(A,F), in the path integral 
representation we have 


P,(A,F)= f 2?[Ft]iP(Ft|Fo)po(ro)d(A[Ft]-A)<5(F,-F) , 

(A6) 

with T*[Ft] = limif_^oo 0^=0 dTi^T/x- In this manner, the 
denominator on the 1. h. s. of Eq. (5), denoted by L for 
convenience, can be rewritten as 


L = J V[Tt] J dFe“Ar)^(r,|ro)po(ro)po(Ft) 

xd(A[F*]-A)d(F,-F) . 


(A7) 


By making use of Eq. (Al) and the property of the delta 
function, we obtain 


L = J V[Tt]ePQ^^d+aArAp(Y^\fo)po{rl) 

xpo(ro)d(A[Ft] - A) . 


(AS) 


The insertion of the delta function aims at constructing 
the joint distribution function of A for the time-reversed 
process and the final phase space point. Particularly, we 
rewrite Eq. (All) as 

L = e^[v[Tt][ dTe^-^^^P{Tt\To)po{ro) , , 

J Jeo (A 12 ) 

xpo(Ft)<5(A[ft]-(-A))(5(f,-F) . 

Recall the path integral representation of the joint dis¬ 
tribution function (A 6 ), Eq. (A12) actually leads to 


L = e^ [ (iFP^(-A,F)po(F^)e®Ar) ^ (A 13 ) 

Jeo 

which is nothing but the numerator on the 1 . h. s. of 
Eq. (5) multiplied by e^. 


Appendix B: Validity of the Dynamic Property (7) 
for the Langevin Dynamics 

We focus on the one-dimensional case here. The gener¬ 
alization to higher dimensions is straightforward, though 
anisotropic effect may emerge. Eor the overdamped 
Langevin dynamics, the generator Ct in the EPE ( 6 ) 
reads [41] 


Ct = -d,{d,Ut + ^-^d,) . (Bl) 

7 

Here 7 is the viscous friction coefficient; Ut = Vt{x) only 
depends on the position, indicating 

4 = A . (B2) 

Notice that dxf{x)£g{x) = - dxg{x)£ f{x) 
for any normalizable functions f{x) and g{x) (so that 
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lim.j:^±oo f{x)g{x) = 0), the transpose operator of Ct 
should be 

Cj = -{-d,Ut + r^d,)d, . (B3) 

7 

Before checking the dynamic property (7), we introduce 
the following two useful relations 

=dj; - iSd^Ut ■ 

=dl - 2pd,Utd, - pdlUt + . 

(B4) 

With these relations in hand, we start to calculate the 1. 
h. s. of Eq. (7) 

=-[dlUt + dM.-P[d,Utf 

(B5) 

+ r^dl - 2dM. - dlUt + ^{d^Utf] 
=U-d.Utd, + r^dl) = Cj . 

7 

Thus, Eq. (7) has been confirmed to be valid for the 
overdamped Langevin dynamics. 

Let’s move on to the underdamped Langevin dynamics. 
The generator Ct reads [41] 


Appendix C: General eqnation of motion for the 
Joint Distribntion Function 

Instead of the definition in the main text (3), lets con¬ 
sider a thermodynamic variable associated with a process 
generally expressed as 

A[rt] = f dtdtWtiTt) + f dtTtdrqtirt) . (Cl) 
Jo Jo 

Here Wt{T) and gt(r) can be arbitrary time-dependent 
functions with respect to the phase space point. We will 
show that if the generator is Ct, the equation of motion 
for Pt{A,r) will be 

dtPt{A,r) = [e- 9 *(r)aA£^g 9 dr) 0 ^ _ dtWt{T)dA]Pt{A,T) . 

(C2) 

It is instructive to first consider a simple case that 
gt(r) = 0. In this case, Pt+dt{A,r) is related to Pt{A,r) 
via the following relation [41] (terms with the magnitude 
of {dt^ are ignored) 

Pt+dt{A, T) =(1 + Ctdt)Pt{A - dtWt{T)dt, r) 

=CtPt{A,T)dt + PtiA - dtWt{T)dt,T) , 

(C3) 


which implies 


Ct = -—dx + dp{da;Ut + 7—) -f 7/3 ^dp , (B6) 

m m ^ 

based on which we can obtain its time-reversed operator 
and its transpose operator 


C\ ——dx + dp{—dxUt + 7 ~) + 7/3 ^dp , 

? . ™ (B7) 

Cj =—dx - {dxUt + l—)dp + 7/3 . 

m m ^ 

Here Ut = f? j2m P Vt{x) depends on both the position 
and the momentum. Again we introduce two useful rela¬ 
tions first 



Then we calculate the 1. h. s. of Eq. (7) 


ePu^de-PU. 

=^{dx-PdxUt) 

m 

+ — + {—dxUt + j—){dp — / 3 —) 

m mm 

+ lP~^dl - 2j—dp - — + 7 / 3 ( —)^ 

^ m m m 

= — dx - dxUtdp - 7—9p -f Jl3~^dp = Cf . 
m m ^ 


(B8) 


(B9) 


Thus, Eq. (7) has also been confirmed to be valid for the 
underdamped Langevin dynamics. 


dtPt{A,T) = [Ct - dtWt(T)dA]Pt{A,r) . (C4) 

This result is familiar to us since it is the precursor of 
the Feynman-Kac formula before performing the integral 
transformation. In order to generalize to the case with 
qt{C) ^ 0, we rewrite the first term in the rightmost of 
Eq. (C3) as 

CtPt{A,r)dt= [ dr'dr"Pt{A,r") 

Jet (C5) 

X 5(r' - r)£((3(r' - r'')dt , 

where £( acts only on Th Now the physical meaning is 
transparent: this is the contribution to Pt+dtiA,T) due 
to motions in the phase space, consisting of both parts 
that come from other phase space points (T" ^ T) and 
leave T (T" = T, this is necessary to correct the second 
term in the rightmost of Eq. (C3), the contribution due 
to the temporal variation of Wt{T)). Based on such an 
interpretation, when a nonzero qtiC) appears, we simply 
modify Eq. (C3) as 


Pt+dtiA,T) = / dT'dV'PtiA - qtiV) + qt(r"),T") 

Jet 

X s{r' - r)£((5(r' - T'')dt + Pt{A- dtwt{T)dt, r), 

(C6) 

owing to the second term in Eq. (Cl). Since the first term 
in Eq. (C6) is already of the order of the magnitude of dt, 
we don’t have to further add modifications like dtWt{C)dt 
or dtqt(r)dt in addition to —gt(r) -|- qt{T''), which will 
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merely result in differences of the order of the magnitude 
of (dt)^. By using the identity e“^/(x) = f(x + a) as 
well as the property of the delta function, we can simplify 
Eq. (C6) as 

+Pt{A-dtWtiT)dt,r) , 

which finally leads to Eq. (C2). 

Now let’s return to the functional in the main text (3). 
By using the identity 


dtdtat{T t) + f 
Jo 

we find that Eq. (3) can be rewritten as 


aT-(r^) - ao(ro) = [ 
Jo 


dtttdratiTt) , 
(C8) 


A[Tt]= [ dtdtat{Tt)+ [ dtTtdriatiTt) -/3Ut{Tt)] . 

Jo Jo 

(C9) 

Comparing Eq. (C9) with Eq. (Cl), we have r(;t(r) = 
at(r) and gt(r) = at(r) — (3Ut{r). Substituting them 
into Eq. (C2), we get the equation of motion in the main 
text (8). 


Appendix D: FT of work (21) for Driven Isolated 
Quantum Systems 


We generally denote the Hamiltonian of an isolated 
quantum system by Jit, whose time-reversal is deter¬ 
mined by T-Lt = , with 0 to be the antiuni¬ 

tary time-reversal operator. To be specific, the quantum 
version of Eq. (21) can be written as 


Y,PAWAn))Poie\n))eP^^ 

n 

=eP'^Y,Pr{-W,\m))po{e-^\m))e^^^ . 

rfi 


(Dl) 


Here |n) (|m)) is an eigenstate of 'Hr {'Hr) with the 
eigenenergy E" (E™); po{Q~^\rh)) (po^^l^^))) is the 
probability that the initial state for the forward (back¬ 
ward) process is measured to be 0“^|m) (0|n)), which 
is obviously an eigenstate of J-Lq {Hq) due to the former 
definitions. Such probability can be evaluated in terms 
of the initial density operator go {go) via 

Po(0”^|to)) =(m|0eo0“^|m) , 
po(0|n)) =(n|0“^po0|«) , 

where [go,'Ho] = [ 00 ,^ 0 ] = 0 due to its structure as¬ 
sumed in the main text. Accordingly, [0po0~^j’Hr] = 
[0 “^Po0i^t] = 0 holds subsequently. Based on the 
two-point projection measurement definition of quantum 


work, the joint distribution functions should be 
Pr{W, |n)) =Y,\{n\Ur,o\m)\^Po{\ni))5{W - + E^) , 

m 

Pr{W, |m)) = ^ \{MUr.,o\n)?Po{\n))5{W - Ef + E^) , 

n 

(D3) 


where po{\m)) (po(|n))) can also be related to go (go) by 
{m\go\m) ((fi|po|h)). 

To prove Eq. (Dl), we again take the characteristic 
function-based approach, which has been widely used in 
the studies of quantum thermodynamics [42, 43]. We 
take the inverse Fourier transformation on both sides of 
Eq. (Dl), and obtain 



dWP^^ Pr{W, \n))po{Q\n))e^^- 

n 


=EE |(n|Wr,o|m)ppo(|m))po(0|«))e^*''+'^^®” 

n m 


n m 

=Tr[e*"«-W.,opoe-*'^^“Wo,r0-ipo0] , 

^ P^(-W, |m))po(0”^|m))e'^^- 


(* + 00 


=EE \{7n\lirfi\n)\'^Po{\n))po{e ^|TO))ed'"+^)^o 

m n 


= ^^(TO|e '■^'^^l(r,ogo\n){n\P'''^°llo,r&goQ ^\m) 

fh n 

=TY[e-^^^'^-Ur,ogoP''^'^Uo,rQQoQ~^] ■ 

(D4) 


Here Ut^t' {Ut,t ') is the time-evolution operator for the for¬ 
ward (backward) process, governed by the Schrodinger 
equation = 1- [ihdtUt^t' = 

= X), I is the identity operator; both p 
and v = p — ifj are generally complex numbers, i.e., p 
and V are unnecessarily their complex conjugates p*, v*, 
and are also independent to each other, owing to the 
arbitrariness of (3. By making use of the algebraic prop¬ 
erties of 0 and the trace (see Ref. [43] for details), espe¬ 
cially Wt,*/ = 0e-*'"’^*0-i = 

Tr[0-ii0] = Tt[A^ and Tt[AB] = Tt[BA] (be careful 
that this may be invalid if A = 0), as well as the com¬ 
mutation relations [po,^o] = 0 and [0 po0~^j^t] = 0, 
we have 

TiiP’''^-Ur.pgoe-^>^'^°Uo,re-^goQ] 

=Ti-[A’''»-^e-^Uo.,rQgoe~"^'^°e~^Ur.pgoQ] 

=Tiie-A-^^''^°Uo.regoQ-A^p'^-Ur,ogoe] 

=TiigoUo,re-^^'^^egoQ-^Ur,oA’'Y 
=Tr[e*"««Wo,r0Po0-'e-*^«-W,,oeo] 
=Tr[e-*'^«-W,,oeoe*""°Wo..0eo0"^] ■ 
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Thus the inverse Fourier transformation of the two sides 
of Eq. (Dl) turns out to be the same. So far, the discrete 
version of the refined FT of work (21) has been confirmed 
to be generally valid for driven isolated quantum systems. 

As an example, if = e“^^”/Zo(/3) and po = 
ZoiP) = Tr[e-'5«o], Zo(/3) = TT[e-^^°] = 
Tr[e-/3«x] = ZriP) (due to Tr[ 0 -M 0 ] = Tr[^t] and 
T-L]. = Hr- This is the quantum analogy of footnote [29]), 
Eq. (D5) will become G{^)/Zr[l3) = G{—fi + i/3)/Z(){f3). 
This means that the quantum CFT holds even for sys¬ 
tems without the time-reversal symmetry, e.g., a charged 
particle subjected to a time-dependent magnetic field. 
This is a generalization of Ref. [43], where [0, H*] = 0 is 
assumed. However, we should be careful that the time- 
reversed Hamiltonian must be , but usually 

not TLr-t- 


Appendix E: Detailed Calculations on the breathing 
Brownian oscillator 

To be consistent with Ref. [33], we use the generating 
function pt(A, x) = Gt{iX, x) instead of the characteristic 
function Gt{fJ-,x). We first present the main result in 
Ref. [33], which focused on an arbitrary process driven 
by the protocol fc* starting from the equilibrium state 
Pq'^{x) = {kol3/2Tr)i By making the Gaussian 
ansatz 


Pt{X,x) 




jt) 

20 A(t) 


(El) 


the Feynman-Kac formula can be self-consistently re¬ 
duced to the following two first-order ordinary differential 
equations 


^A(i) = —, (E2) 


<^A(t) 


7 P7 


3Afct 
2 il^xit) 


(E3) 


with the initial condition to be '!/ 1 ’a( 0 ) = 1 and (()a( 0 ) = 
(/3A:o)“^. To perform perturbative analysis, it is con¬ 
venient to define gx{t) = I3kt(l)x{t)/tpx{t). Based on 
Eqs. (E2) and (E3), it can be checked that gx{t) satisfies 
the following Raccati equation 


Accordingly,ln'!/)A(T) can be expressed as follows up to 
the first-order accuracy 


1 / / ^ A jx X r 

InV’A('r) =-^In—- —(1 --) / 
2p kq 4p p Jq 


+ C>{p) . 

(E5) 

In particular, for the protocol fcj = feg + Ave have 

Now let’s start to calculate the generating function of 
the work distribution for the nonequilibrium initial state 
Pq{x) = {k^P' /2TT)^e~^d'kox'^ ^ Using Eq. (46), we obtain 


(e 


-AW\ 


=g-/3AF 


ipoG 

^+oo 


/ -l-oo 

dxpriP — A, x) 

-oo 


Pojx) 


/3'fco[#-A(T)]3 /•+“ + 


2TTPkr4>p-xiT) 7 - 


/. 




= + "V'/3-A(r) 


(E7) 


where all the quantities with an overline must be asso¬ 
ciated with the time-reversed protocol fc* = k^-t- After 
the first-order approximation, Eq. (E7) becomes 

=\ ;§)4-aW + 

= ln'0^^^(r) - ^(1 - j,)9^plx^T) + 0{Y) , 

(E 8 ) 


where = — 7 KA/( 2 fcQ/ 3 ), and the fluctuation- 

dissipation relation ln'i/’A^^('r) = —PAF + 
is used. Since the generating function in the form 

ln(e“'^^) = —X{W) + must correspond to a 

Gaussian distribution centred at (VF) and with variance 
cr^, we finally obtain the results (47) in the main text. 


Appendix F: Derivation of the Sagawa-Ueda 
equalities (53) and (54) from the Unified IFT 


<?a(0 = - —[<?A(t) - 1] + ^gx[t)[l - ^gx{t)] , (E4) 

7 kt p 

with the initial condition to be gxY) = 1- Once gx{t) 
is determined, 'ipxi'r) = {e~^^) can be obtained via 
ln'0A(T) = -{X/P) dtgx{t)kt/2kt. If we regard fc* 
in Eq. (E4) as a small quantity, the zeroth-order so¬ 
lution will simply read g^P(t) = I, while the first- 
order correction should be gP (t) = 7(1 — X/P)kt/2p. 


We consider the extended phase space as the direct 
(Cartesian) product of the phase spaces of the system 
and the measurement device. The phase space point in 
such extended space can be denoted by Et = (Ft, y), with 
Ft and y to be the components of the system and the 
device respectively [44] . The unified IFT of the composite 
system reads 


Pri^r) 

Po(Eo) 


)= 1 - 


(FI) 
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This is the key relation that we will use. 

To derive Eq. (53), we choose 

-/3C/o(ro) 

p„(E„)=p(!/|r„) . 

By substituting Eq. (F2) into Eq. (FI), we obtain 


(F2) 


as long as the DB condition holds. To distinguish two 
contributions to A.Stot' (i) the correlation between the 
system and the device, and (ii) the entropy production 
of the system and in the medium, we define the mutual 
information at the initial and the final stages 


/. = In , ^/ = In ■ 

Po{^o)p{y) pr{Tr)p{y) 


(F7) 


,•^0^ /3[Q+C/«(rO-C/o(ro)]-ln[p(y|ro)/p(y)]\ 

^zy{p) ^ 


1. (F3) 


Then we can obtain another decomposition identity of 
Astot as follows 


Here the DB condition, and thus Asm = PQ has been 
assumed. To simplify Eq. (F3), we make use of the 
first law kFpt] = Q[rt] + f/r(r,-) — ?7o(ro) as well 
as the expression of the free energy difference AF^ = 
—P~^ \a.[Zy{P)/ZQ{P)], and further define the state func¬ 
tion of the mutual information = ln[p(?/|ro)/p( 2 /)]. 

Combining these relations we obtain 

{e-^(w-AF)-i^ = 1 . (F4) 

To derive Eq. (54), we choose p“(Er) to be the real 
final distribution priZ^-). So we get 

= 1 , (F5) 

which is nothing but the EPI for the composite system. 
Concretely, we have 

Astot =/3Q - lnpr(S.r) + Inpo(So) , (F6) 


Astot = PQ - lnp.r(rr) + Inpo(ro) - If + h = a - AI , 

(F8) 

which leads to 

(e-<^+A/) ^ . (F9) 

It should be mentioned that po(S) and p“(S) are tac¬ 
itly assumed to be nonzero for any E to guarantee the 
validity of the two Sagawa-Ueda equalities in this ap¬ 
pendix. However, this assumption is usually not satis¬ 
fied for error-free measurements or a rigorously localized 
initial state, thus the two Sagawa-Ueda equalities may 
break down in these cases [11]. Recent researches have 
shown that by adding an extra modification term on the 
exponential, the Sagawa-Ueda equalities can be general¬ 
ized to be applicable to the feedback control processes 
with error-free measurements [45]. 
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